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Abstract 

Wc argue a new classification scheme of long gamma-ray bursts (LGRBs) using the morphology of the 
cumulative light curve of the prompt emission. Wc parametrize the morphology by the absolute deviation 
from their constant luminosity (ADCL) and derive the value for 36 LGRBs which have spectropic redshifts, 
spectral parameters determined by the Band model, 1-second peak fluxes, fluences, and 64-msec resolution 
light curves whose peak counts are 10 times larger than background fluctuations. Then we devide the sample 
according to the value of ADCL into two groups {ADCL < 0.17 and ADCL > 0.17) and, for each group, 
derive the spectral peak energy E p - peak luminosity L p correlation and the Fundamental Plane of LGRBs, 
Q^,) ■ which is a correlation between the spectral peak energy E p , the luminosity time Tl (= Ei SO /L p where E- lso 

is isotropic energy) and the peak luminosity L p . We find that both of the correlations for both groups 
are statistically more significant compared with ones derived from all samples. The Fundamental Planes 
with small and large ADCL are given by L p = 10 52 - 53±0 01 (-E P /10 2 - n keV) 1 - 84±0 03 (T L /10 - 86 sec) a29±0 08 
with xl = 10.93/14 and L p = 10 52 - 98±0 - 08 ( J Ep/10 2 - n keV) 1 - 82±0 - 09 (T L /10 - 86 sec) - 85±0 - 27 with xl = 7.58/8, 
respectively. This fact implies the existence of subclasses of LGRBs characterized by the value of ADCL. 
Also there is a hint for the existence of the intermediate- ADCL class which deviates from both fundamental 
planes. Both relations arc so tight that our result provides a new accurate distance measurement scheme 
up to the high rcdshift universe. 
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Introduction also imply the existence of subclasses and if we can iden- 
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■ tify them it would not only reduce the ambiguity of the 
In spite of the discovery of the observational diversity correlation but also give some information on the origin 

\^ . of the afterglow, parent galaxy, and environment proper- of LGRBs. 

ties of gamma-ray bursts (GRBs) in Swift era, the prompt In this paper, we show the statistical significance of 

emission studies over the last 10 years have not succeeded the correlation improves by considering the morphology 

in revealing telltale GRB subclasses beyond the short- of the cumulative light curve of the prompt emission, 

hard and long-soft classification (Kouvcliotou ct al. 1993). Specifically, first, we define the absolute deviation from 

■ Except for some unusually low luminous events which are their constant luminosity (ADCL) of the cumulative light 
j_h " often associated with Type Ic Supernovae (e.g., 980425, curve. Then we derive the correlation both for small and 
5^ ■ 060218), all long gamma-ray bursts (LGRBs) seem to have large ADCL groups and they are shown to be tighter than 

similar properties and no subclasses. one derived by all LGRBs. We also find that dimmer out- 

In Tsutsui et al. (2009), we found a correlation between hers of Tsutsui et al. (2011) have intermediate ADCL 

the spectral peak energy E p , the luminosity time Xl and values, 
the peak luminosity L p (Fundamental Plane 1 of LGRBs). 

However, the correlation has a relatively large number of 2. Data Description 
outliers (~ 30%) and they deteriorate the usefulness of 

LGRBs as standard candles which probe the expansion Since the discovery of a correlation between the peak 

history of universe (Tsutsui et al. 2011). The outliers may energy E p and the isotropic energy Ei so (Amati et al. 



2002), several correlations between E p and other bright- 

The Fundamental Plane is a two-dimensional family first discov- / v , n , i nm , rii • i j . i onni \ 

j . , , . . . . . , . . , , ness (Yonetoku et al. 2004; Ghirlanda et al. 2004a), were 

crcd between velocity dispersion, mean surface brightness and v . 

luminosity of elliptical galaxy (Djorgovski & Davis 1987). The proposed. There are also some studies to improve these 
term is used for various two-dimensional families of other ob- correlations by adding third parameters of GRBs (Liang & 
jects, such as galaxy clusters (Schaeffer et al. 1993), black hole Zhang 2005; Firmani et al. 2006; Tsutsui et al. 2009; Yu 

activity (Merloni et al. 2003) and X-ray gas in normal galaxies et al 2009). These correlations are thought to be im- 
(Diehl & Statler 2005). , , / , , • , ,i , r • 

portant not only to investigate the nature of their emis- 
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sion mechanism but also to measure distance up to the 
high redshift universe where no other distance indica- 
tors have been observed (Ghirlanda ct al. 2004b; Schaefer 
2007; Kodama ct al. 2008; Liang ct al. 2008; Cardone et al. 
2009; Tsutsui et al. 2009). 

In spite of their attraction, there have been many debate 
against these correlations (Nakar & Piran 2005; Band & 
Preece 2005; Butler et al. 2007; Rossi et al. 2008; Collazzi 
& Schaefer 2008; Shahmoradi & Nemiroff 2011). 

Most of these discussions, however, are based on incom- 
plete data without redshifts, and data with many system- 
atic errors which conceal the nature of GRBs. Therefore, 
we do not take these works into account in this paper. 

We have investigated origins of systematic errors of 
the spectral - brightness correlations (Yonetoku et al. 
2010; Tsutsui et al. 2010), and succeeded to remove them 
(Tsutsui et al. 2011). As a result, we have found that 
the Ep-T^-Lp correlation possibly divides LGRBs into 
at least three subclasses. In this paper, we study a new 
classification scheme based on a relationship between the 
residual from the Ep-T^-Lp correlation and long term 
pulse shape within the whole duration of bursts. 

Before we study the classification of LGRBs, here we 
update and recompile the database of Yonetoku ct al. 
(2010) and Tsutsui et al. (2011). There are more than 
120 GRBs with known redshifts and spectral parameters 
of time- integrated spectrum up to October 2011. However 
more than half of their spectral parameters determined by 
not the Band model (Band et al. 1993), but by the Cut-off 
power law model (Pendleton et al. 1997) because of the 
lack of number of high energy photon and/or the limited 
energy range of detectors, which tends to overestimate 
spectral peak energies (Kaneko et al. 2006; Krimm et al. 
2009; Tsutsui ct al. 2010). Using these biased data, the 
difference between subclasses is smoothed out with their 
systematic errors (Tsutsui et al. 2011). Therefore we do 
not deal with these data from the beginning in this pa- 
per. Then, we use only 44 LGRBs with known redshifts, 
spectral parameters determined by the Band model, ob- 
server frame 1-second peak fluxes, fluences, and 64 msec 
resolution light curves. 

For the GRBs detected by GBM/ Fermi, the spectral 
analysis is performed with the software package RMFIT 2 
(version 3.3rc8) and the GBM Response Matrices vl.8. 
We analyze the time-integrated spectrum using CSPEC 
data, from 8 keV to 40 MeV with 1.024 sec temporal res- 
olution, following the guidance of the RMFIT tutorial 3 . 

For other GRBs, there are spectral parameters of indi- 
vidual GRBs reported by different observation teams and 
different authors, so that we choose the data fitted by the 
Band model with the smallest uncertainties. Because we 
did not know the systematic difference between E p esti- 
mated with the Band model and E p with the CPL model 
in our previous studies, the data in this paper is not ex- 
actly the same as the data in Yonetoku et al. (2010) and 
Tsutsui ct al. (2011). 



The KONUS team often reports not 1 sec peak fluxes 
but 64, 128, 256 msec fluxes, so that we convert these 
peak fluxes to 1 sec peak fluxes with 64 msec light curves 
obtained from the archive of all the GCN/KONUS Notices 
and Light Curves 4 (Yonetoku et al. 2010). 

For all data reported with errors at 90% confidence 
level, we convert it to errors at 68% one by multiplying 
1/1.645. In table 6, we summarize redshifts, spectral pa- 
rameters, observed 1-second fluxes and fluences with de- 
tector's energy band and their references of 44 LGRBs. 

From these parameters, bolometric 1-second luminosi- 
ties (£p b f s ) and isotropic equivalent energies (i?i so ) be- 
tween 1-10,000 kcV in GRB rest frame are calculated by 
following equations: 



Lpl = 47T<iL-Pp i obs 



rlO,000/(l + z) 
Jl/(l+z) 



E x N(E)dE 



Eu 



J^N(E)dE 
(erg s _1 ), 



l + z 



Sobs * 



Jz™ x Ex N(E)dE 
(erg). 



(1) 



(2) 



where N(E) is the Band model in units of 
photons cm~ 2 s -1 keV -1 (Band et al. 1993) and c?l 
is the luminosity distance in units of cm. We should note 
that ip b i's are calculated within 1-second in observer 
frame, which must be converted to rest r-second peak 
luminosities Lj^'s (Yonetoku et al. 2010; Tsutsui ct al. 
2011). In this paper, we adopt r = 2.752 second for all 
GRBs. 3). 

There is another parameter characterizing a time scale 
of GRBs. Firmani et al. (2006) used the high signal time 
Tqa5 to improve the E p -L p correlation. Here we adopt 
another time scale, the luminosity time Tl, first used by 
Tsutsui ct al. (2009) defined as, 



E[ S 



(3) 



where an error of Tl is calculated from the error propa- 
gation formula without the cross-term. 

In this paper, we analyze light curves of LGRBs to iden- 
tify subclasses of LGRBs. For secure analysis of light 
curves we use only the bursts whose background fluctu- 
ation - peak counts ratio are less than 0.1. This reduces 
the number of GRBs from 44 to 36. Then wc analyze 36 
bursts in following sections. In table 7, wc summarize the 
intrinsic properties of the 36 GRBs. 

3. Fundamental Plane of Long Gamma-Ray 
Bursts 

The time scales which improve the E p -L p correlation 
have been studied in the past. Some obtained positive re- 
sults (Firmani et al. 2006; Tsutsui et al. 2009) and the oth- 
ers negative results (Rossi et al. 2008; Collazzi & Schaefer 



2 http:/ /fcrmi. gsfc.nasa.gov/ssc/data/analysis/ 

3 http:/ /fcrmi. gsfc.nasa.gov/ssc/data/analysis/user/vc_rmfittutorial.'pdf http:/ /gen. gsfc.nasa.gov/konus_grbs. html 
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2008). There are several reasons that the previous stud- 
ies could not establish the tight correlation between the 
spectral peak energy, the time scale, and the peak lu- 
minosity : systematic errors of peak energies estimated 
by the CPL model, systematic errors of peak luminosi- 
ties estimated in 1 seconds of observer, and contamina- 
tions of short GRBs with extended emission or other out- 
liers. They have all disturbed the estimation of the best 
fit function and their dispersion so that such tight corre- 
lations have never been reproduced without the studies 
of the systematics as mentioned above (Yonetoku et al. 
2010; Tsutsui et al. 2010; Tsutsui et al. 2011). 

Here we briefly describe the method developed in 
Tsutsui et al. (2011). We assume a linear correlation be- 
tween E p , Tl, and L p in a logarithmic form as: 



logL p (E p ,T h ) = A + Blog^ 
+C\og 



10 2n keV 



10 086 sec 



(4) 



where A, B, and C are free parameters of the model. If 
we adopt C = 0, and C = — 1, this is identical with the 
Ep-Lp and E p —Ei so correlations, respectively. 
The corresponding \ 2 function is given by : 



(5) 



logL^-logLpOEp^TL,;) 



The factor 2C in front of <J\ og L comes from the fact that 
the definition of Tl includes L p . 

Because the \ 2 value strongly depends on the existence 
of outliers which deviate from the Gaussian distribution, 
we first optimize the Lorentzian merit function, given by 



M{A,B,C,t)=J2HI + 1^1 



(7) 



with <Ti n t = 0. Because the Lorentzian merit function is less 
affected by outliers which do not follow Gaussian distri- 
bution, we use the parameters optimizing the Lorentzian 
merit as the tentative values to eliminate outliers and to 
estimate robust a- m t value. 

To estimate robust dint; we use the robust standard 
deviation given by, 

_ median [|logLJ$ - logL p (E Pti ,T w )\] 

,JRSD = 0^745 ' (8) 

and then we obtain robust a m t as follows : 



'int — U RSD 

N 



-l^{(l + 2C)< ipi+ B 2 < 



E p ,i +^ 2 °fogT L . 1 } 



(9) 



Now that we have tentative sets of parameters and <j- u 



we can compute, for each sample, t = \z%\ and the two- 
tailed P-valuc from the t distribution with (TV — 3) degrees 
of freedom. We adopt a threshold probability as Pth = 
Q/N for all GRBs, and the GRBs which have P-value 
smaller than P tn are regarded as outliers and eliminated 
from the following chi square analysis. Here Q is arbitrary 
number less than 1. If it is not mentioned we always use 
Q = 0.2 in this paper. Therefore we now have data sets 
without outliers and dint , and we can perform the ordinary 
chi square analysis to estimate the best fit parameters and 
their uncertainties. 

First we apply this method to our updated sample, 
both the whole sample and the platinum sample with 
aE p /E p < 0.1. In Tsutsui et al. (2011), using only the plat- 
inum sample we obtained the very tight correlation with 
six outliers. The reason for only the platinum sample is 
to distinguish possible outliers or subclasses from the ob- 
servational errors. As a result in Tsutsui et al. (2011) we 
suggested that there might be three subclasses in LGRBs. 
In figure 1, we plot the E p -Tl—L p diagram obtained by the 
present platinum sample. The platinum events consistent 
with the Ep-T^-Lp correlation, that is, on plane events, 
are marked with red circles, and outliers with green tri- 
angles. In table 1, we summarize the best fit parameters, 
fraction of outliers (N out /N), intrinsic dispersion (<7i n t) 
and reduced chi square (xt) both for the whole sample and 
the platinum sample. As table 1 shows, the best fit values 
of parameters, <7i n t, and the fraction of outliers (N ou t/N) 
depend on the sample. We see that a lnt for the whole 
sample is 2.6 times larger than that for the platinum sam- 
ple. While the outliers exist for both samples so that the 
distributions are not Gaussian for both samples. This fact 
suggests that strong non-gaussianity, perhaps bimodality 
of the scatters around the best fit function exists. As fig- 
ure 1 indicates, there might be a subclass of GRBs (green 
triangles in the upper left corner) which might be another 
correlation nearly parallel to the Ep-T^—Lp correlation. 
This situation is similar to PoPI and PoPII Cepheid vari- 
able correlations that are parallel due to the difference of 
metallicity. In practice, however, it is sometime difficult 
to remove all the outliers as shown in Monte Carlo simula- 
tions in Appendix A in Tsutsui et al. (2011) even if we use 
the outlier detection technique based on robust statistics 
. This might be the reason for the best fit parameters and 
dint for platinum sample in table 1 being different from 
the ones in Tsutsui et al. (2011) although the tightness 
of the correlation is unchanged. This fact suggests that 
there is a large fraction of outliers close to the best fit 
function. If we can remove these outliers using other ob- 
servational properties in advance, the Fundamental Plane 
of LGRBs will be more statistically significant and ro- 
bust. The purpose of this paper is to find what is this 
other observational property like the metallicity in PoPI 
and PoPII Cepheid variable correlations. 

In the following sections, we introduce a new parame- 
ter of LGRBs, the absolute deviation from their constant 
luminosity (ADCL). We show that dividing LGRBs ac- 
cording to ADCL- value significantly improves both of the 
Ep-Lp correlation and the Fundamental Plane of GRBs. 
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That is, we show that ADCL is a new observational prop- 
erty what we have been seeking for. 

Table 1. Fitting results for our new sample. 



sample 



all 



platinum 



best fit 



(52.63,1.63,-0.08) 



(52.52,1.86,0.31) 



CTint 



0.13 



0.05 



N out /N 



8/36 



11/23 



by ourselves. The most popular duration of GRBs is Tqq 
defined by Kouvcliotou et al. (1993). Tgo is the time dur- 
ing which the cumulative counts increase from 5% to 95%. 
However, we find that Tgo tends to lose weak long tails of 
bursts which are essential to classify the bursts. Thus, we 
44.se Tgg to define ADCL instead of Tgg. 



20.9/26 



3.81/10 



M IS 
£ o 









* mr 




J* 


m 






f s 









\0~ 2 10 _1 10° 

(E p /10 2 - 71 [keV]) 1 - 86 (T L /10' 



10 1 1(T 
fl.86 [sec]) 0.31 



Fig. 1. The E p — Tl~L p diagram of updated platinum sam- 
ple. The green triangles represent outliers. 



4. Absolute Deviation from Constant Luminosity 
(ADCL) 

To quantify the absolute deviation from their constant 
luminosity (ADCL), we analyze 64 msec resolution light 
curves of BATSE, KONUS, Swift, and GBM. We summa- 
rize the energy band of each detector in table 2. 

Table 2. The energy band of light curves for each detector 



BATSE 


KONUS 


Swift 


GBM 


20-1000 kcV 


50-200 kcV 


15-150 kcV 


8-1000 kcV 



We determine the background intervals before and af- 
ter bursts for BATSE, Swift, and GBM light curves and 
make linear fit of background using least square method. 
For KONUS light curves, there are no pubic data before 
bursts so that we use the background intervals only after 
bursts. Then we subtract the backgrounds from 64 msec 
light curves interpolating or extrapolating the fit across 
the whole duration of bursts. From these background sub- 
tracted light curves, we determine the duration of bursts 
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Fig. 2. An example of the cumulative light curve of GRB 
990506. The cumulative light curve is plotted with a red 
dashed line, and the ordinary light curve with a gray solid 
line. The vertical lines show the start and end time of Tgg. 
The interval within Tgg is indicated by an arrow. 

Because the observed distribution of GRBs spans the 
wide range of redshift, from 0.0085 to 8.2 up to present, 
and the cosmological time dilation changes the time scales 
of bursts by factor (1 + z), the measurement of temporal 
properties of bursts must be robust to the change of time 
resolution of light curves. Besides, recent results from hy- 
drodynamical simulations of jet imply that the long term 
variability of GRBs come from the activity of the central 
engine, while the short term variability from the propa- 
gation of jet (Morsony et al. 2010) and we should decom- 
pose long and short term variability in some way. For 
these reasons, we use cumulative light curves to quantify 
the difference of bursts which are much simpler than ordi- 
nary light curves (McBreen et al. 2002; Varga et al. 2005). 
Figure 2 shows an example of the cumulative light curve 
of GRB 990506. The cumulative light curve normalized 
with total counts is plotted with a red dashed line, and 
the ordinary light curve normalized with peak counts with 
a gray solid line, respectively. The short scale structure 
in the ordinary light curve is canceled out in the cumula- 
tive light curve, and then the analysis of cumulative light 
curves is much robust to a cosmological time dilation and 
choice of the time resolution. 

To compare pulse shapes of bursts with various dura- 
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tions, we normalize the time scale of bursts with Tgg and 
then all bursts have normalized cumulative counts (7?°™ 
from 0.01 to 0.99, and normalized time tf° rm from to 1. 
Figure 3 shows normalized cumulative light curves of the 
platinum samples in table 7. On-plane events are plotted 
with red lines, outliers with green lines in figure 3. The 
black dotted line is the line of the virtual source with the 
constant luminosity. The brighter and dimmer outliers 
deviate from the constant luminosity line, while many of 
on-plane events cluster around the constant luminosity 
line. 
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Fig. 3. The normalized cumulative light curves for the plat- 
inum samples. On-plane events are plotted with red solid 
lines, outliers with green solid lines. The black dotted line is 
the line of the virtual source with the constant luminosity. 

To quantify these difference, we introduce a new param- 
eter of GRBs : the absolute deviation from their constant 
luminosity (ADCL). We define the absolute deviation 
from their constant luminosity : 



ADCL = 



\Cf 



0.01-0.98 x f5 



AT bi] 



(10) 



where the number of the bin (iVbi n ) is different from burst 
to burst. Figure 4 shows the distribution of ADCL. The 
distribution of ADCL in figure 4 indicates that there is a 

gap around log(ADOL) 0.75 (ADCL ~ 0.17). From 

this result, we divide the LGRBs into two subclasses : one 
is small- ADCL events with ADCL < 0.17 and the other is 
long tailed ones with ADCL > 0.17. Figure 5 shows exam- 
ples of cumulative light curves of small- ADCL (left) and 
long tailed (right) events, respectively. Although it is diffi- 
cult to quantify the statistical significance of the existence 
of the gap, we will show that the gap becomes more clear 
when we study the relationship between the ADCL and 



-1.2 -1.0 -0.8 -0.6 
Log(ADCL) 



-0.4 



Fig. 4. The distribution of ADCL. 
log(ADCL) 0.75. 



There is a gap around 



residual from the Fundamental Plane of LGRBs. This is 
very similar to the discovery of the short - hard and long - 
soft classification by Kouvcliotou ct al. (1993). Only from 
Tgo, the long and short classification is not so clear that 
they added the spectral hardness. Then they could iden- 
tify the bimodal distribution of the long-soft and short- 
hard GRBs clearer. In our case, the residual from the Ep- 
TL-Lp correlation corresponds to the hardness. In Figure 
8. we can identify two suclasses (blue squares and red 
circles) in the ADCL-residual plane. 

In table 7, we summarize the intrinsic properties and 
ADCL of the 36 GRBs. 

4-1- Type I Fundamental plane 

In Tsutsui et al. (2011) and §-3, we have found that 
there are probably three classes of LGRBs : on-plane 
events, brighter outliers, and dimmer outliers. In section 
4, we introduced a new parameter ADCL and found that 
the most of on-plane events have smaller ADCL, while 
most of the brighter and dimmer outliers have relatively 
larger ADCL. In this section, we derive the best fit func- 
tion and dispersion of the E^-T^-Lp correlation for small- 
ADCL and large- ADCL GRBs separately, using the same 
method in section 3. 

We used only the platinum data with ue v ID v < 0.10 in 
Tsutsui et al. (2011) and section 3 to avoid the contam- 
ination of outliers and subclasses because there are no 
prior knowledge about outliers and subclasses of LGRBs. 
Now ADCL is a candidate for another observational prop- 
erty to distinguish outliers and subclass of LGRBs. Using 
ADCL we might eliminate brighter and dimmer outliers in 
advance. Therefore we do not need to restrict the sample 
to the platinum one from now on. The platinum sam- 
ple was needed only to find ADCL so that we will use 



G 



R. Tsutsui et al. 



the whole sample in the following analysis although larger 
observational errors compared with the platinum sample 
exist. 

First, we apply the outlier elimination technique to the 
small- ADCL GRBs with ADCL < 0.17 and find the D p - 
T L -L p correlation of small-ADCL GRBs becomes tightest 
around r = 2.752 seconds. The best fit relation is given 
by, 

1.84±0.03 



10 



o2.o3±0.01 



L„ 



10 2n keV 



10 oj 



0.29±0.08 



(11) 



with xl = 10.93/14 and a- m t = 0. The power-law index of 
Tl in equation. (11) is different from that of the equation. 
(12) of Tsutsui et al., (2011). This is because there were 
some contamination of long tailed GRBs for the lack of 
population of LGRBs. Hereafter we refer to the equation 
(11) as the Type I Fundamental Plane. 

In figure 6, we plot the data of the small- ADCL events 
with ADCL < 0.17 and the best-fit function of the Type 
I Fundamental Plane. The small-ADCL events used to 
derive the best-fit function were marked with red circles 
and the events regarded as outliers with green triangles. 
The solid line represents the best-fit function of the Type 
I Fundamental Plane. 
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Fig. 6. The Type I Fundamental Plane of small-ADCL 
events (ADCL < 0.17). The red circles represent the small- 
-ADCL event used to derive the Type I Fundamental Plane 
and the green triangles represent outliers eliminated from chi 
square analysis. The solid line represents the equation (11). 



4-2. Type II Fundamental Plane 

We find that largc-ADCL GRBs form the other 
Fundamental Plane which is almost parallel to Type 



I Fundamental Plane but slightly brighter than it. 
Furthermore, we find that all GRBs which have large 
Tl > 12 sec deviate from other events, and we eliminate 
them from following analysis. Then we apply the out- 
lier elimination technique to the long tailed events with 
T L < 12 sec and estimate the best fit parameters and res- 
olution time scale of long tailed events. We find the E v ~ 
T^-Lp correlation of the long tailed events also becomes 
tightest around 2.7 seconds. Adopting 2.752 seconds as 
the resolution time scales of large- ADCL events, we ob- 
tain 



10 



52.98±0.080 



L„ 



10 2n keV 

T L 
10 086 sec 



1.82±0.093 



0.85±0.26 



(12) 



with xl = 7.58/8 and Omt = for long tailed GRBs. 
Hereafter we refer to the equation (12) as the Type II 
Fundamental Plane. The Type II Fundamental Plane is 
about 2.5 times brighter than the Type I Fundamental 
Plane. 

In figure 7, we plot the data of large- ADCL events 
and the best-fit function of the Type II Fundamental 
Plane. The events used to derive the best-fit function 
were marked with blue squares and the events regarded 
as outliers with green triangles. The solid line represents 
the best-fit function of the Type II Fundamental Plane. 

Figure 8 shows the relationship between residual 
from the Type I fundamental plane and ADCL of all 
events. Each symbol represents same as figure 6 and 
7. Dimmer outliers apparently have intermediate ADCL 
value. Figure 8 might imply that there are some other 
classes of LGRBs. 

In section 3 we showed that the D p -Tl-L p correla- 
tion for whole sample has significant intrinsic dispersion 
(cint = 0.13 for whole sample while 0.05 for platinum sam- 
ple, see Table. 1). If we divide the whole sample into two 
subclasses from their ADCL values, a- m t = as shown in 
this chapter. Thus we can conclude that the E p —Tl,—L p 
correlation is certainly improved by introducing ADCL. 

5. Discussion 

In this section, we would like to show the validity of 
the E p -Tl-L p correlation. First one might ask that Tl 
might introduce spurious correlation because Tl includes 
L p in its definition. We checked this by assuming the re- 
lation as logLp(L p ,L; iso ) =A'+B' log(L p /10 2 n kcV) + 
C' log(L iso /10 53 - 40 sec). We determined A', B and 
C using the same data and the same method. If there is 
no spurious correlation, we should obtain B = B / (1 + C) 
and C = C/(l + C) where A, B and C are defined in 
Eq. (4). We obtained that B = 1.42, C' = 0.23 with 
xl = 10.93/14. B = 1.84, C = 0.29 so that B/(l +C) = 
1.42, C/(l + C) = 0.22 for small-ADCL sample and that 
B = 0.98, d = 0.46 with xl = 7.59/9. B = 1.82, C = 0.85 
so that D/(l + C)=0.98, C/(l + C) =0.46 for large ADCL 
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Fig. 8. The relationship between the residual from Type I 
fundamental plane and ADCL of all events. Each symbol 
represents same as figure 6 and 7. 

sample. These values show that the spurious correlation 
is very small even if it exists. Therefore which one should 
be used is a matter of taste. 



Next we compare the E p -T^- 



L p correlation with E p ~ 



j p correlation for each subsample with the same statistical 
We first derive the E p -L p correla- 



samplc. As table 1 and 3 indicate the E p -±j p 
has the larger systematic dispersion than the Ep—ij^-^p 



method in section 3. 
tion with the whole sample as discussed in section 3. In 
table 3 we summarize the best fit parameters, fraction of 
outliers (N out /N), intrinsic dispersion (a- lnt ) and reduced 
chi square (xt) both for the whole sample and platinum 

L„ correlation 
T L -L F 

correlation when we use the whole sample, the whole sam- 
ple. 

Secondly we divide the GRBs into two groups by 
ADCL, and derive the E p -L p correlation for each sample. 
Here we should note without the E p -T^-L p correlation 
it was not possible to find ADCL and to recognize the 
existence of the subclasses. However now we know that 
two subclasses exist so that we check the E p -L p correla- 
tion for each subclass. In table 4 and 5, wc summarize 
the best fit parameters, fraction of outliers (N out /N), in- 
trinsic dispersion (dint) and reduced chi square (xt) f° r 
each correlation and each subclass. As table 4 and 5 in- 
dicate, the difference in parameters of the E p -L p correla- 
tion between small- ADCL events and large- ADCL (long 
tailed) events is smaller than that of the E p ~Tl~L p corre- 
lation. Therefore we can not distinguish them even if we 
use our outlier rejection technique since the difference of 
the amplitude A is almost the same (~ 0.2 in logarithm 



E p -L p correlation while 



Therefore E p -T L - 



-0.45 for E p -T L - 



L p correlation) 



two subclasses. In case of small- ADCL events, the E p -L p 
correlation has slightly larger dispersion than the E p -T^- 
L p correlation, although the E p -L p correlation has one 
more outlier. Taking into these facts, we can insist that 
the E p -T\,-L p correlation is surely more statistically sig- 
nificant than E p -L p correlation for small- ADCL events. 



In case of long tailed events, however, it is difficult to say 
which one is more significant. The E p -L p correlation has 
two more outliers, but smaller dispersion than the E p — 
T^-Lp correlation. To show the validity of the E p -Tl-L p 
correlation for long tailed events, we need larger number 
of data. 



Table 3. Fitting results for the E p -L p correlation with the 
whole sample. 



sample 


best fit 


Cint 


N ou t/N 


xi 


all 


(52.64,1.62) 


0.18 


6/36 


27.6/28 


platinum 


(52.57,1.71) 


0.15 


6/23 


12.9/15 



Lp correlation is needed to distinguish 



6. Summary 

In this paper we defined ADCL which characterizes 
the cumulative light curve of the prompt emission of long 
GRBs. We divided the events into two groups accord- 
ing to the value of ADCL. Then we derived Fundamental 
Planes for the two groups separately and it was found that 
they are tighter and statistically more significant than one 
derived from the whole sample. Although wc introduced 
the new parameter ADCL to reduce the dispersion of the 
Fundamental Plane, we found dividing the events into two 
groups also improves the E p -L p correlation. This fact 
supports the validity of introducing ADCL. These tighter 
relations would be helpful to use them as standard candles 
to probe cosmological expansion history. Also our anal- 
ysis may imply the existence of subclasses of long GRBs 
classified by ADCL. 

The discovery of these possible subclasses of GRBs re- 
minds us of the discovery of the two separate subclasses of 
Cepheids by Baade. After Baade's discovery, the accuracy 
of the distance measurement by the Period-Luminosity 
relation of Cepheids was drastically improved. We can 
expect that the distance measurement by the Ep-T^-Lp 
relation of GRBs will also be improved, and we are study- 
ing the distance measurement by the Type I Fundamental 
Plane of GRBs in near future. 
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Table 4. Fitting results for small- ADCL events. 



correlation 


best fit 


N out /N 




xi 




(52.54,1.81) 


5/21 





15.55/14 


Ep-T^-L p 


(52.53,1,84,0.29) 


4/21 





10.93/14 


Table 5. Fitting results for long tailed events. 


correlation 


best fit 


AU/N 


Cint 


xl 


Ep-Lp 


(52.73,1.78) 


6/15 





4.672/7 


Ep-T^-Lp 


(52.98,1,82,0.85) 


4/15 





7.58/8 
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Fig. 5. left : An example of a cumulative light curve of small-ADCL GRB. right : An example of a cumulative light curve of long 
tailed GRB. 
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Fig. 7. left : The Fundamental Plane of large- ADCL events. The blue squares represent the events with Tl < 12 sec used to estimate 
the Type II Fundamental Plane and the green triangles outliers eliminated from chi square analysis. The solid line represents the 
equation (12). right : The relationship between residual from the equation (12) and Tl of the long tailed events. Each symbols 
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Table 6. Spectropic redshifts, spectral parameters determined by the Band model, observed fluxes and fluences with energy band 
of detectors, and references of 44 LGRBs 
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" Taken from the J. Greiner's GRB table and references therein (http://www.mpe.mpg.de/ jcg/grbgen.html). 
b Observed peak photon (or energy) fluxes between E\ to E2 keV. 
c Observed fluences between E3 to E4 keV. 

d Detectors that provide spectral parameters and light curves (SAX=BeppoSAX, K=Konus, BAT=BATSE, HE= HETE-II, S= Swift, 

SK = Swift and Konus, SW= Swift and WAM) 
e References for spectral parameters, peak fluxes and fluences : (1) Amati et al. (2002) ; (2) Yonetoku et al. (2004) ; (3) Sakamoto et al. 

(2005) ; (4) Krimm et al. (2006) ; (5) Golenetskii et al. (2006) ; (6) Golenetskii et al. (2007) ; (7) Krimm et al. (2009) ; (8) Golenetskii 

et al. (2008a) ; (9) Golenetskii et al. (2008b) ; (10) This work ; (11) Golenetskii et al. (2008c) ; (12) Golenetskii et al. (2009) ; (13) 

Golenetskii et al. (2011a) ; (14) Golenetskii et al. (2011b) ; (15) Golenetskii et al. (2011c) 
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Table 7. Intrinsic properties of 36 LGRBs 
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a I : Small- ADCL GRBs consistent with the Type I Fundamental Plane. 

II : Long tailed GRBs consistent with the Type II Fundamental Plane. 

O : Outliers. 
6 E p = (l + z)E° ha 

c Integrated between 1-10,000 keV in GRB rest frame. The flat ACDM universe with Hq = 70 km s _1 Mpc and S^m = 0.30 is assumed. 



